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An investigation of the effects of the incidence angle on the aerothermodynamic 
environments of the Mars Science Laboratory has been conducted. Flight conditions 
of peak heating, peak deceleration and chute deploy are selected and the effects of 
the angle of attack on the aerodynamics and aerothermodynamics are analyzed. The 
investigation found that static aerodynamics are well behaved within the considered 
range of incidence angles. Leeside laminar and turbulent computed heating rates 
decrease with incidence, despite the increase in the leeside running length. 
Stagnation point was found to stay on the conical flank at all angles of attack, and 
this is linked to the rapid flow expansion around the shoulder. Hypersonic lift to 
drag ratio is limited by the heating rates in the region of the windside shoulder. The 
effects of the high angle of incidence on the dynamic aero at low Mach remains to be 
determined. Influence of the angle of attack on the smooth-wall transition 
parameter indicates, that higher angle of attack flight may result in delayed 
turbulence onset, however, a coupled analysis, involving flight trajectory simulation 
is necessary. 


Nomenclature 

L/D = lift to drag ratio 

M e = boundary-layer edge Mach number 

R - radius, m 

Re = Reynolds number 

U,v = velocity, m/s 

T = temperature, K 

a = angle of attack, degrees 

p = density, kg/m 3 

Subscripts 
b = base 

e = boundary-layer edge 

max - maximum over the surface 

n - nose 

oo = free- stream 

0 = boundary-layer momentum thickness 
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I. Background 

The present analysis is a segment of a trade study to determine the range of acceptable lift to drag ratio 
for Mars Science Laboratory (MSL) entry. MSL is a lifting Viking- derived vehicle with a heatshield cone 
half-angle of 70 degrees (Figure 1). The range of acceptable incidence angles is bound by the severity of 
the aeroheating environment at hypersonic speeds, and the aerodynamic stability at low speeds and at 
parachute deploy. Since MSL, due to its large size, will encounter turbulent heating during its decent, the 
magnitude of the turbulent heat rates as a function of the angle of attack is evaluated. The location of the 
stagnation point and the peak laminar heating is determined for angles of attack between 10 and 26 degrees. 
It is found that for most of the conditions, considered in this analysis, the stagnation region remains on the 
cone flank and does not migrate to the shoulder of the vehicle. This is especially important during the high 
heating segment of the entry, as the movement of the stagnation point to the shoulder may result in the 
augmentation of the shoulder heating beyond the thermal protection system (TPS) capability. The physical 
phenomenon that occurs at high 
hypersonic conditions and limits the 
location of the stagnation region to 
the cone flank even for the incidence 
angles greater then 20 degrees is 
found to be a consequence of the 
sonic line location on blunt cones 
with small shoulder radius. 

The present analysis indicates that 
predicted laminar and turbulent 
heating over most of the surface of 
the vehicle decrease with the increase 
in incidence within the considered 
range. The heat rate at the windside 
shoulder increases, as could be 
expected, and, for some conditions 
exceeds the heating over the rest of 
the vehicle. 

The peak predicted value of the 
smooth wall transition parameter 
Re 0 /M e over the surface is found to 
decrease with the increase in 
incidence. Based upon the completed Figure 1. MSL entry vehicle geometry 

CFD it may be possible to develop a 
relation between the length of the 

turbulent heat pulse and the trim angle of attack. While the present data indicates that the incidence 
increase may reduce the duration of the turbulent heating, a higher fidelity analysis is required. Analysis of 
the effect of the angle of incidence on the length of the turbulent heating pulse would require a flight 
simulation to properly construct the new trajectory, since the higher L/D vehicle will generally fly along a 
different trajectory, then that of a lower L/D vehicle. 

The principle objective of this analysis was to determine the maximum acceptable L/D for MSL. One of 
the concerns is the behavior of the static aero at incidence angles beyond 16 degrees. As the present 
investigation indicates, the static aero at the peak heating, peak deceleration and parachute deploy flight 
conditions are very nearly linear between 10 and 26 degrees incidence and do not limit the selection of L/D 
for the mission. 



2 


II. Aerothermal environments 


A. Computational Method 

Aerothermal environments were calculated using LAURA program 1,2 . This program is a parallel multi- 
block code that is extensively used in entry vehicle calculations. The code can solve Euler, full Navier- 
Stokes and Thin Layer Navier-Stokes flow- fields using an upwind point- and line-implicit relaxation. In the 
present analysis the solutions of the reacting Thin Layer Navier-Stokes equations are employed. Fully 
turbulent solutions were also calculated using Baldwin-Lomax algebraic turbulence model. The gas (eight 
specie model: C02,C0,C,N2,N,N0,N2,02,0) is assumed to be in thermo-chemical non-equilibrium. The 
wall catalysis is enforced by requiring recombination of C02 at the vehicle’s surface. This condition is 
presently used in the aeroheating analysis to correlate the test data. The heating rates are calculated at 
radiative equilibrium wall temperature. 

Two conditions, selected for these calculations represent the peak heating/peak deceleration region of 
the 0517 peak heat rate MSL trajectory. At the peak heating point, which occurs at around Mach 23.4 
flight, the laminar solutions at 10, 18 and 26 
degrees incidence illustrate general trends. A 
more detailed analysis of the influence of 
the angle of attack on hypersonic 
aerothermodynamics is developed at Mach 
16 flight condition. At Mach 16 laminar 
solutions are calculated between 1 0 and 26 
degree incidence in 2 degree increments. 

Turbulent solutions are also computed at 
several angles of attack. 

In determining the aerodynamic heating 
of Viking-derived shapes it is common to 
calculate the forebody flow fields separately 
because the flow over the shoulder is 
supersonic for all typical combinations of 
incidence and the free stream Mach for 
which aeroheating is important. This makes 
it possible to decouple the forebody flow 
from the wake-flow. Since the boundary 
layer edge Mach number on the wind-side 
shoulder decreases with the increase in 
alpha, it is important to verify that the 
forebody decoupling is permissible for the 
range of alpha considered in this analysis. 

B. Influence of Incidence on Computed 
Heatrates 

Figure 2 and Figure 3 show the variation 
of boundary layer edge Mach number in the 
vehicle’s pitch plane. Wind-side and leeside 
shoulders of the vehicle are on the left and 
right sides of the plots respectively. The 
figures demonstrate that the shoulder flow is 
supersonic for all incidence angles at Mach 
23.4 and 16. Forebody flow is, therefore, 
free from the influence of the wake at these 
conditions, and can be calculated separately. 



Figure 2. Pitch plane M e (Moo= 16) 



x (m) 


Figure 3. Pitch plane M e (Moo= 23.4) 
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The influence of the vehicle incidence angle on the laminar heat rates and their distributions is shown m 
the Figures 4 and 5. The wind-side shoulder heating increases with alpha, while the heating over most of 
the rest of the forebody decreases. Interestingly, similar trend exists in the turbulent heating distribution 



Figure 4. Mach 23.4 LAURA laminar heat rates 



Figure 5. Mach 16 LAURA laminar heat rates 



Figure 6. Mach 16 LAURA turbulent heat rates 


(Figure 6). Turbulent heating levels over most of 
the body decrease with increasing incidence, 
despite the increase in the leeside surface 
streamline running length. Solutions at higher 
incidence angles indicate that the wind- side 
shoulder heating may exceed leeside flank heating. 
The leeside flank heating rates presently determine 
TPS sizing, and this is acceptable for moderate 
angles of attack. 

C. Effect of the Incidence on Surface Pressures 

Variation of the surface pressure in the plane of 
symmetry as a function of the angle of attack is 
shown in the Figure 7. Increase in incidence angle 
moves the center of pressure away from the 
centerline of the vehicle, which may result in lower 
compressive forebody deformations. It is expected 
that MSL will experience static aeroelastic 
deformations, which, depending on their scale, may 
attribute to changes in the aerodynamics. Unlike 
most of its predecessors, which had flown non- 
lifting trajectories at incidence angles near zero, 
MSL is an actively guided lifting vehicle, and the 
trim characteristics are an important part of the its 
aerodynamic performance. Due to the precision 
landing requirement 3 , MSL is particularly sensitive 
to factors, such as the aeroshell deformations, 
which introduce aerodynamic uncertainties during 
flight. While the present approach is to design the 
forebody such as to minimize deformations 
altogether, it may be advantageous to investigate 
high incidence operation as a possible way to 
reduce deformations for the future missions of this 
class. 



Figure 7. Mach 16 pitch plane surface pressures 
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D. Stagnation point movement 

It was believed that the increase in incidence angle past the cone half angle will result in the stagnation 
point movement to the shoulder of the vehicle and aggravation of the vehicle’s aerothermal environment. 
We find that, at hypersonic regimes the stagnation point stays on the cone, does not transition to the 
shoulder, and its movement along the cone originator is very nearly linear with respect to the angle of 
attack as its varied over the range, considered herein. The variation of the wall shearing stress in the 
vehicle’s plane of symmetry (Figures 8, 9 and 10) illustrates the motion of the stagnation region as the 
incidence angle is varied. 

Due to the small shoulder radius, windside shoulder flow is mostly supersonic for all cases, indicating 
that the flow stagnation will occur on the cone flank for a wide range of incidence angles. While the shear 
stress distribution may suggest asymptotic behavior this is not so. The flow will stagnate on the shoulder 
once the incidence angle, sufficient to move the sonic line further off the cone is reached, but the incidence 
angle required for this to occur is outside the interrogated range (Figure 11). 




Figure 8. Pitch plane shear, Mach 23.4 


Figure 9. Pitch plane shear, Mach 16 



x(m) 



Figure 10. Pitch plane shear, Mach 2 


Figure 11. Influence of incidence angle on the sonic line, Mach 23.4 
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III. Influence of angle of attack on turbulent transition 

Because of the aeroshell size and the ballistic coefficient, MSL is going to experience turbulent heating 
during entry 4 . An experimental program, involving LaRC Mach 6 Air tunnel, CUBRC LENS and GALCIT 
T5 facilities has been under way to bound turbulent augmentation and to determine the onset of transition 5 . 
While MSL is being designed for a fully turbulent heat pulse, better understanding of the length of the 
actual turbulent heat pulse and the effects of the angle of attack on turbulent heating is desirable. 

To determine the onset of turbulence and the likely length of the turbulent heat pulse a set of laminar 
solutions along the trajectory is computed, as outlined in the Table 1. Angle of attack for these solutions 
was constant and equal 13 degrees. Laminar solutions yield the peak values over the surface of the 
transition parameter at each condition. It should be pointed out, that the values of the maximum Ree/M e 
over the surface are approximate in the sense that they are estimated from the CFD. A number of numerical 
and user uncertainties feed into this, such as the calculation of the momentum thickness and boundary layer 
thickness by the code, and determining the proper location of the maximum Re 0 /M e by the analyst. 

Table 1. Conditions along 0517 peak heat rate trajectory 


Time 

(sec) 

Altitude 

(m) 

(m/s) 

Poo 

(kg/m 3 ) 

Tco 

(K) 

668 

5.126E+04 

5.521E+03 

3.004E-05 

137.3 

686 

3.412E+04 

5.242E+03 

2.671E-04 

160.4 

695 

2.692E+04 

4.766E+03 

6.833E-04 

170.4 

703 

2.212E+04 

4.089E+03 

1.139E-03 

175.3 

717 

1.833E+04 

2.861E+03 

1.968E-03 

180.6 

727 

1.727E+04 

2.230E+03 

2.216E-03 

181.9 

746 

1.601E+04 

1.545E+03 

2.229E-03 

182.5 

840 

9.263E+03 

5.504E+02 

4.026E-03 

189.8 


A number of researchers believe that the turbulent transition occurs when the parameter Re 0 /M e reaches 
300 somewhere on the surface of the vehicle. Other researchers believe that a number closer to 200 is more 
appropriate. We have consistently seen numbers between 350 and 550 in the tests, conducted to date 5 . The 
solutions, completed in this study provide some insight into the behavior of the smooth wall transition 
parameter Re 0 /M e at various angles of attack. Figures 13 and 14 indicate decrease of the peak Re 0 /M e with 
increase in angle of attack. The only deviation is the 10 degree incidence Mach 23.4 case, for which the 
peak Re 0 /M e doesn’t occur at the centerline of the vehicle (Figure 15) unlike 18 and 26 degree incidence 
cases for the same Mach (Figures 16, 17). The ten degree incidence case is the only case where we see this 
behavior, and it has to do with a different flow structure, with a stagnation point on the spherical cap. This 
case aside, it is evident that an increase in incidence tends to reduce Re 0 /M e over most of the vehicle. The 
possible consequences are the reduction of the peak turbulent heating and the delay of turbulent transition 
until later in the trajectory. Note that the discontinuous behavior of the parameter in the vicinity of the 
stagnation point (Figures 10, 1 1 and 14) has to do with the calculation of the {0/0} type indeterminacy. 

It is important to note, however, that the present analysis is only a first order look at the phenomena. 
Typically, a higher L/D vehicle will fly a different trajectory from that of a lower L/D vehicle, hence, a 
more detailed study would be needed to determine the connection between the angle of attack and the 
transition onset. Analysis of the effect of the angle of incidence on the length of the turbulent heating pulse 
would require a flight simulation to properly construct the new trajectory. 
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Figure 12. Peak Re 0 /M e during entry 
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Figure 13. Re 0 /M e in plane of symmetry (Mach 23.4) 



Figure 14. Re 0 /M e in plane of symmetry (Mach 16) 


Figure 15. Mach 23.4, 10 degree incidence 



Figure 16. Mach 23.4, 18 degree incidence 
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IV. Static Aerodynamics 

Hypersonic aerodynamics of Viking-derived entry vehicles are dominated by the forebody. Even at 
high angles of attack, when the flow on the windside of the backshell is partly attached (Figure 18) the 
predicted contribution of the aftbody is not significant (Table 2). The analysis of the effects of the 
incidence on the hypersonic 
aerodynamics is, therefore, carried out 
on the forebody only. Solutions at Mach 
16 and Mach 23.4 flight conditions 
indicate linearity of the aerodynamics 
over wide range of incidence angles 
(Figures 18-23). 

Fow supersonic aerodynamics, on 
the other hand, could be expected to be 
influenced by the aftshell to a greater 
degree, since the aftshell contribution is 
a function of the type of 1/M 2 (see, for 
example, 6) presently applied to axial 
force coefficient only. Recirculation 
region changes in size and location on 
the aftshell, as the angle of attack is 
varied. The aerodynamics of the aftshell 
are complex and are quite non-linear in 
incidence. Computed M=2 aero for the 
forebody only and the entire aeroshell is 
shown in Figures 21 and 22. We find 

that the predicted aftbody contribution to the pitching moment and F/D is still rather small, regardless of 
the incidence angle. 



Figure 18. MSL wake, Mach 22.6, 15 degree incidence 




Figure 19. Hypersonic L/D Figure 20. Hypersonic pitching moment 

(Mach 16, Xcg/D - radial CG offset) 


Table 2. Mach 22.6 aero at 15 degrees incidence 



CA 

Cm 

CN 

Aeroshell 

1.46 

-0.0437 

0.0499 

Heatshield 

1.46 

-0.0436 

0.0491 

Aftshell 

O(E-05) 

O(E-05) 

O(E-04) 
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Figure 21. Supersonic L/D 


Figure 22. Supersonic pitching moment 


V. Summary and Conclusions 

The influence of the angle of attack on the aero thermodynamic environment of the Mars Science 
Laboratory was investigated as part of the study to bound the acceptable lift to drag ratio. Angle of attack 
was found to reduce both the computed laminar and turbulent heating rates over most of the heatshield 
while increasing the windside shoulder heating. It is found that the stagnation region stays on the conical 
flank despite the increase in incidence angle past the cone half-angle. This is due to the rapid flow 
expansion around the shoulder, which effectively limits the sonic line from moving further downstream. 

No non-linearities were observed in the flow field structure over the interrogated range of the free- 
stream conditions. Peak value of the smooth wall turbulent transition parameter Re 0 /M e was found to 
decrease with increase in incidence, despite the increase in the lee-side surface streamline running length. 

The static aerodynamics are linear at both the supersonic and hypersonic Mach over the entire range of 
incidence angles. The computed aftshell contribution to the total aero is small, and the influence of the 
incidence angle on the wake doesn’t significantly affect the static aero. 
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